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ABSTRACT 



^0 ' in circumstellar disks, which are commonly referred to as "transitional" or "cold" disks, and it is the first scenario to check before 



Context. The presence of um'esolved binaries on sub-arsecond scales could explain the existence of optically thin inner holes or gaps 



making any other assumptions. 

Aims. We aim at detecting the presence of companions inside the inner hole/gap region of a sample of five well known transitional 
disks using spatially-resolved imaging in the near-lR with the VLT/NAC0/S13 camera, which probes projected distances from the 
Q . primary of typically 0.1 to 7 arcsec. The sample includes the stars DoAr21, HD 135344B (SAO 206462), HR4796A, TCha, and 

in ' TWHya, spanning ages of less than 1 to 10 Myr, spectral types of AO to K7, and hole/gap outer radii of 4 to 100 AU. 

^ ' Methods. In order to enhance the contrast and to avoid saturation at the core of the point-spread function (PSF), we use narrow-band 

C^ , filters at 1 .75 and 2.12 fim. The "locally optimized combination of images" (LOCI) algorithm is applied for an optimal speckle noise 

removal and PSF subtraction, providing an increase of 0.5-1.5 mag in contrast over the classic method. 

Results. With the proviso that we could have missed companions owing to unfavorable projections, the VLT/NACO observations rule 
out the presence of unresolved companions down to an inner radius of about 0'.' 1 from the primary in all five transitional disks and 
^ , with a detection limit of 2 to 5 mag in contrast. In the disk outer regions the detection limits typically reach 8 to 9 mag in contrast and 

^~7 ■ 4.7 mag for T Cha. Hence, the NACO images resolve part of the inner hole/gap region of all disks with the exception of TW Hya, for 

\l ' which the inner hole is only 4 AU. The 5cr sensitivity profiles, together with a selected evolutionary model, allow to discard stellar 

companions within the inner hole/gap region of T Cha, and down to the .sM&.?te/Zar regime for HD 135344B andHR4796A. DoAr21 is 
v^ I the only object from the sample of five disks for which the NACO images are sensitive enough for a detection of objects less massive 

r~^. . than ~ 13 Mji,p that is, potential giant planets or low-mass brown dwarfs at radii larger than ~ 76 AU (0'.'63). 

f^ ' Conclusions. These new VLT/NACO observations further constrain the origin of the inner opacity cavities to be owing to closer 

or lower-mass companions or other mechanisms such as giant planet formation, efficient grain growth, and photoevaporation (for 
DoAr21andHR4796A). 



Key words, protoplanetary disks - planetary systems - circumstellar matter - stars: individual: TCha, HD 135344B, DoAr21, 
HR 4796A, TW Hya - instrumentation: adaptive optics - stars: imaging 

1. Introduction ISkmtskie et al.l [l990l) identified several sources with spectral 

energy distributions (SEDs) that showed a deficit of mid-infrared 

The transition from the gas/dust-rich actively accreting optically fl^^ and a rise into the far-IR. These were interpreted as a sign 

thick disks surrounding T-Tauri and Herbig Ae/Be stars to the ^f just clearing, which is expected from disk evolution as the 

gas/dust-poor optically thin debris disks is currently not well un- just grows and settles or is removed with the gas by viscous ac- 

derstood. Of the several processes proposed for disk evolution ^^^^^^^ ^^ ^jher dispersal mechanisms such as photoevaporation 

and dissipation, the formation of planets remains the most exci- ^nd tidal truncation because of close companions, both stellar 

ting (e.g., | Meveretal. || 2007i and references therein). Early stu- ^nd substellar. More recent surveys of young stellar clusters with 

dies of young stars dating back to IRAS (e.g., [ Strom et al. | 1989|: gp^^^^ revealed a small population of these sources in several 

near by star-forming regions (d < 400 pc) based on p hotom- 

Send offprint requests to: S. Vicente, svicente@rssd.esa.int etry dSicilia-Aauilar et al. 2006; Muzerol le et al.l 2010h. spec- 



Based on observations collected with the NAOS-CONICA (NACO) troscopy (.Brown et al. .2007: .Furlaii et alj 120091: lOliveira et al.l 



instrument at the VLT in Paranal Observatory, operated by the European |2010h or both JSicilia-AVuilar et al.l l2008nMerin et al.l |2010|) . 
Southern Observatory, under programme ID 079.C-0914(A) 
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Table 1. Observing log of the 2.12 jum VLT/NACO images of the five stars with transitional disks and corresponding PSF-calibrators 



Star 
&PSF 


Date 


Grade' 


DIT 

(s) 


NDIT 


Exp.time- 
(min.) 


Seeing-' 
farcsec) 


Airmass 
(at start) 


(ms) 


StrehP 

(%) 


FWHM<5 
(mas) 


V' 
(mag) 


Notes'* 


DoAr21 
CD-24 13079 


2007-06-18 


B 


10.0 
4.0 


4 
1,3 


8.0 
0.73 


0.68 
0.75 


1.530 
1.547 


2.55 
2.14 


9.3 

15.2 


129.1 
85.0 


14.01±0.03 
9.86 


sat. 
not sat. 


HD 135344B 
HD 136961 


2007-05-09 


B 


1.0 
1.0 


12 
10 


2.4 
1.67 


1.75 
1.85 


1.062 
1.071 


2.73 
2.14 


13.5 
14.5 


82.5 
84.6 


8.708±0.017 
6.754 




HD 135344B 
HD 144156 


2007-06-20 


B 


2.0 
1.0 


6 
10 


2.4 
1.67 


0.77 
0.79 


1.051 
1.076 


1.07 
1.23 


25.1 
26.1 


74.4 
72.3 


8.708±0.017 
8.50 


close comp. 


HR4796A 
HD 109536 


2007-06-14 


B 


1.0 
1.0 


10 
1 


2.0 
0.17 


NA 
NA 


1.069 
1.092 


1.67 
1.88 


35.9 
27.0 


69.2 
71.9 


5.78 
5.127 




TCha 
CD-78 512 


2007-05-23 


B 


25.0 
12.0,2.0 


2 
1,6 


10.0 
1.0,0.8 


0.85 
1.16, 1.89 


1.735 
1.739, 1.740 


1.37 
1.28,0.76 


4.8 
7.8,7.6 


116.2 
97.5,97.5 


11.86 
10.14 


ellip. 
ellip. 


TWHya 
HD 94889 


2007-05-12 


A 


5.0 
1.0 


5 
5 


5.0 
0.83 


1.46 
1.10 


1.033 
1.053 


1.27 
0.59 


19.0 
23.9 


78.4 
75.7 


11.07 
8.80 


close comp. 


TWHya 
HD 101636 


2007-05-23 


B 


5.0 
0.3454 


5 
5 


5.0 
0.29 


0.98 
0.86 


1.020 
1.026 


0.94 
0.09 


18.3 
30.8 


77.9 
71.0 


11.07 
9.41 





' the grade of each observing block (OB) is set by the on-site VLT observer based on the user-requested ambient conditions (seeing, sky transparency, airmass) vs. the actual observing 

conditions. A grade 'A" means that the conditions were fully within specifications, "B" that they were mostly within specifications, while "C" means they were out of specifications 

and the OB should be repeated if possible. 
2 on source = DIT X NDIT X number of OBJECT frames 
^ average astronomical site monitor seeing at start in the visible (0.5 fim) and at the zenith (airmass- 1). This value can difli'er significantly from the NACO image quality, which is 

usually better. "NA" means that no seeing data was available for this night. 
'' average coherence time of the atmosphere at 0.5 /im calculated by the real time computer (RTC) on real data. Short values of tq (< 6 ms) represent a seeing not stable in time 
^ computed with eclipse in the final processed image 
^ of the star in the processed image adopting a Moffat fitdng to the radial profile. The full width at half-maximum (FWHM) of the difl:"racdon-limited PSF for the VLT at 2. 1 2 /jm is 

65 mas 

' optica l magnitude re trieved from SIMBA D: 1) DoAr21 - ICieza et al.l 12001 , 2) HD 135344B - IHeg et alj t2000h , 3) HR4796A - IBarrado Y Navascuej|200g) , 4) TCha - ITanner et al] 

<2007|) , 5) TW Hya - ITorres et alj 120061 ) 
^ DoAr21: PSF saturated in the central pixels (4 to 7 pixel range); TCha: elliptical along x-axis; PSF calibrators HD 144156 and HD 94889 have close companions at d=l'.'33, PA=20° 

and d=2'.'2, PA=260°, respectively. 



The y are currently refer red to as 'transitional' or 'cold' disks 
after lBrown et al.l (l2007h . because of their lack of emission from 
warm dust. The reported fraction of transitional disks can vary 
from a few to 50% depending on definition and the method of 
classification. A consensus in nomenclature is still lacking and, 
as demonstrated in Men'n et al. (2010), a significant fraction of 
transitional disks classified with photometry turns out to be not 
transitional when using spectroscopy. The general definition of 
a transitional disk is a young disk (< 10 Myr) with an opti- 
cally thin inner region (the "opacity hole/gap") surrounded by 
an optically thick outer disk (beyond the hole/gap outer radius, 
Rhoie/gap.out)- The mid-infrared emission originates in the inner 
edge or "w all" of the truncate d disk that is directly illuminated 
by the star (ICalvet et alJl2005l) . Gapped disks that still retain an 
inner component of warm optically thick dust, typically of less 
than a few AU in width, ar e also referred to as 'pre-transitional' 
after lEspaillat et al] (l2007h . Most transitional disks are still ac- 
creting, so that gas is being transported through the inner cleared 
hole/gap onto the star. In some sources a small amount of micron 
or sub-micron dust coexists with the gas in the hole/gap region 
and causes an excess over photospheric fluxes in the near-IR. 

The detection of tidal gaps in some disks has confirmed 
the inference of the inner holes from unresolved data. These 
gaps and holes have already been resolved with high-resolution 
contrast imaging and interferometry in the optical, infrared 
and at (s ub)millimeter waveleng t hs in a handful of objects: 
LkCal5 dThalmann et alJ 1201 ft iPietu et all l2006l). T WHy a 



(iHughesetalJ I2007h. LkHg 330 (iBrown et all 120081 l2009l) . 
GMAur ('Dutrev et al.' '2008; ' Hughes et al.1 l2009h . HR4796A 
(Schneider etal. 1999, 2009), SR21N, WSB60, DoAr44, 



eling of transitional disks suggests inner cavi ties in the dust 
distributi on of 3 to 100 AU in extens i on (e .g., D'Alessio et alJ 
2005.: .Naiita et al J 120071: [Brown et al.l 120071: ICieza et al.1 120 1 Ol: 
Olofsson et al.ll201lT) . Although several mechanisms have been 



HD13 5344B (lAndrews et aP 120091: iBrown et alJ l2009l) . and 
TCha dOlofsson et al.ll20111) . The direct images and SEDs mod- 



proposed to explain the origin of inner disk holes or cleared gaps, 
the existence of previously unresolved close binaries has been 
confirmed in se veral transitional disks dlreland & KrausI [20081 : 
iPott et alll201Cll) and is the first scenario to check before making 
any other assumptions. 

This paper presents near-infrared high-spatial resolution 
narrow-band imaging of five well known transitional disks - 
DoAr21, HD135344B (SAO 206462), HR4796A, TCha, and 
TWHya - and is aimed at detecting close companions within 
their inner hole/gap region. The VLT/NACO observations and 
data reduction are presented in § |2] The results, including the 
PSF-subtracted images and the sensitivity curves are shown in 
§ |3] where we also refer to the mass limits and orbit constraints 
of possible companions. In §|4]we tentatively classify the five 
disks on basis of the mechanisms that potentially create their in- 
ner opacity hole/gap. We also analyze the results for each target 
independently and discuss how the constraints derived from the 
NACO images complement or support previous studies of com- 
panions that made use of other techniques. The conclusions of 
this study are given in § |5] 

2. Observations and data reduction 

2.1. Observations 

The target stars were observed with NACO, an adaptive op- 
tics supported high-resolution imager and s pectrograph installed 
on ESO's Very Lar ge Telescope in Chile ( iRousset et al.l 120031 : 
iLenzen et al.ll2003h . The observations were performed in ser- 
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Fig. 1. LOCI optimization grid for a full width at half-maximum 
(FWHM) of 6 pixels overlaid on the 1 .75 yum PSF-subtracted im- 
age of HD 135344B (09-05-2007). The frame is 512 x 512 data 
pixels (13.27 mas/pixel) and the diameter of the inner circle is 
~ 52 pixels. The FWHM input used in the LOCI pipeline was 
adjusted to the average FWHM of the data for each target (for 
most of the data FWHM=6 pixels and for some FWHM=9 and 
10 pixels). The relevant input parameters for LOCI are not the 
exact annulus sizes but the surface area of the ring subsections A 
(see explanation in section I2.2l i. which in our case corresponds 
to 70 times the PSF area. The algorithm adjusts the ring sizes 
to match this area condition (70 PFS cores/segment) and conse- 
quently the ring sizes change slightly for different FWHMs and 
for different radial positions (because the farther out, the more 
segments). 



vice mode in several runs between May and June 2007 with two 
narrow-band filters, NB.2.12 (2.122 + 0.011 ^^m) and NB_1.75 
(1.748 + 0.013 jum), to avoid saturation in the core of the PSF. 
The images were collected with the highest spatial resolution 
camera S13 (13.27 mas/pixel, 14" x 14" FoV) and with the 
NAOS visible wavefront-sensor and the visible dichroic. The 
camera 'Double RdRsRd' readout mode and the 'HighDynamic' 
mode (15 000 ADU Full Well Depth) were selected to observe 
these bright stars. The PSF standard stars were chosen to be close 
in spectral type and parallactic angle to the target stars. They 
were observed immediately after the science targets to match 
sky conditions as well as possible. More details on the observa- 
tions can be found in Table [1] The observing conditions of the 
different runs were mostly poor (bad and highly variable seeing, 
clouds, winds), resulting in a low adaptive optics (AO) perfor- 
mance. The analyzed datasets barely met the required specifica- 
tions; part of the data had to be discarded. Only one dataset was 
fully within specifications. Hence, future datasets obtained un- 
der more favorable atmospheric conditions will certainly further 
improve the results. 

2.2. Data reduction and the LOCI algorithm 

We adopt a strategy different from the "classic" PSF subtraction 
approach that takes advantage of the stability of the quasi-static 
speckle pattern of the AO-corrected PSF. Each individual image 



(DIT) of the target is dark-subtracted and flat-fielded using stan- 
dard reduction procedures with IRAF and the NACO pipeline, 
which is based on the eclipse library (Devillard 1997). For each 
of them, the stellar PSF speckles are attenuated by subtracting an 
optimized PSF obtained using the "locally opti mized combina- 
tion o f images" (LOCI) algorithm detailed in iLafreniere et al.l 
(l2007h . Conceptually, this data processing technique uses the 
fact that high-contrast imaging PSFs from a given instrument - 
even when taken with different observational setups - share com- 
mon speckle structures. If an observing strategy is arranged in 
a way that a set of images with similar speckle patterns can be 
taken while a faint companion is present in one image and absent 
in others, one can build a 'best estimate' of the PSF for a given 
image from the set of reference images (PSF images). This 'best 
estimate' is subtracted from the target frame and a significant 
gain in sensitivity can be achieved. 

This approach was initially conceived in the context of angu- 
lar differential imaging or ADI (see demonstration in Figure 8 of 
Lafreniere et al.lf2007l and its comparison with the iMarois et al.l 



20061 ADI results), where a long sequence of observations is 
taken with the Cassegrain rotator turned 'off'. Through the se- 
quence, a companion at a given field position will rotate through 
the field as the parallactic angle changes but the average speckle 
pattern caused by telescope and instrument aberrations typically 
remains stable. With this sequence, a high-precision PSF frame 
can be built for each image frame from all other frames taken at 
sufficiently different angular positions (to avoid self-subtraction 
of a companion). This technique can be extended to any set of 
images where a putative companion is absent or significantly 
suppressecQ. That means it can also be applied to a conventional 
data set (field rotation compensated for by the Cassegrain rotator 
following) using PSF stars that are separately recorded as we do 
in this paper. 

The heart of the LOCI algorithm can be described quite easi- 
ly: a specific target frame is divided into "optimization subsec- 
tions" and the best linear combination of reference images is 
found for every optimization subsection to represent the corres- 
ponding target frame section. The idea to split the target image 
into subsections takes into account that the correlation between 
target and reference (PSF) images typically varies with its field 
position. A crucial step in the LOCI algorithm is the definition 
of a proper set of reference images. In angular differential ima- 
ging data for e.g., the displacement caused by field rotation has 
to be large enough to avoid self-cancellation. This complication 
does not exist in our case. We simply used another set of PSF 
star observations with an identical instrumental setup. But we 
exploit the fact that NACO PSFs show a strong degree of mirror 
symmetry and also add the mirror-flipped versions of each PSF 
frame to the set of reference frameo For the NACO data we 
adjusted the surface area of the optimization subsections A (as 
seen in the grid in Figure [TJ to contain approximately 70 "PSF 
co res", i.e., we setLafren iere's parameter A^^i to 70 in Equation 1 
of ILafreniere et al.l (|2007|) : 



iFWHM\- 
NattI- 



(1) 



2 / 

where the FWHM is the full width at half-maximum of the PSF 
and Na the number of "PSF cores" that fit in the optimization 



E.g. by spectral differential imaging (ADI). See in lArtigau et alJI2QQ8l) the application 
of LOCI to the Gemini Near-Infrared Coronographic Imager (NICI) data, a dual-channel 
imager exploiting the methane absorption feature for planet search. 

" Note that obviously LOCI provides a higher gain in sensitivity for ADI data because 
there is no additional slew to acquire the PSF star and furthermore the pupil does not rotate. 
Therefore, there is a better match of the quasi-static speckle pattern for the reference and 
target frames. 
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Fig. 2. High-contrast residual images after PSF removal with the LOCI algorithm for both filters NB_1 .75 and NB_2. 12. Each panel 
is a3'.'3x3'.'3 sub-image centered on the target star The bottom set of panels is a zoom-in of the inner l"xl" region for the NB_2. 12 
filter showing the PSF residuals at small distances from the primary. No speckles were found to be achromatic (i.e., no evidence for 
a real object) and present in the two epochs of HD 135344B (~1.5 month apart) and TWHya (just a few days apart) for which only 
the best dataset is displayed. North is up and east to the left. 



subsection. We did not pursue an extensive optimization for this 
parameter but kept it large enough for convenient calculation 
times and to be on the safe side concerning self-cancelation. 
Figure[T]shows the annulus sizes and geometry of the optimiza- 
tion subsections A for a FWHM of 6 pixels. 

After subtracting the reconstructed PSF from every frame, 
the images were registered and median-combined. The remai- 
ning speckles add up incoherently in the final image, while any 
companion would add up coherently. The final images or PSF- 
subtraction residuals are shown in Figure |2l 



(achromatic) or diff'erent datasets when available, as in the case 
of TWHya (just a few days difference) and HD 135344B (-1.5 
month). A careful visual analysis of the images discarded any 
achromatic speckles that could point to a tentative detection of 
nearby sources. The best images of HD 135344B (20 June 2007) 
and TWHya (12 May 2007) are shown. We note that even 
though our data are not ideal for LOCI application (poor AO per- 
formance, low contrast, and speckles are less dominant), the re- 
sults are still significantly better than those obtained with "clas- 
sic" PSF subtraction (see Figure[3]l 



3. Results 

Figure |2] shows the final images after the data reduction and 
analysis with the LOCI high-contrast technique described above. 
Each tile (3'.'3 x 3'.'3) is centered on the target star and each row 
displays the residual images in one of the two narrow-band fil- 
ters. The bottom set of panels is a zoom-in of the inner 1" x 1" 
region for the NB_2. 12 filter The optimal suppression of speckle 
noise and the improvement in the sensitivity limit allows us to 
search for eventual companions. Some speckles remain in the 
images after LOCI and are caused by the non-optimal observa- 
tional strategy. Electronic ghosts aligned with the X and Y axis 
and row saturation effects typical of NACO were present in many 
of the processed frames of the target stars and PSF-calibrators. 
These residuals can be discarded by blinking between the 1 .75 
and 2.12 fim images, or between different datasets if they were 
acquired close in time and hence are not affected by the or- 
bital motion of a possible companion. A real companion looks 
like a speckle that does not move between the different filters 



3.1. Detection and mass limits for putative stellar/substellar 
companions 

In order to quantify the upper limits on mass of potential com- 
panions, we compute the differential magnitude profile as a func- 
tion of the separation from the primary star. They are obtained 
by calculating the 5cr RMS over concentric annuli (1 FWHM 
width) of the intensity on the residual maps (PSF-subtracted) 
displayed in Figure |2] These resulting sensitivity curves are 
shown for each filter and image in Figure |3] The black and pink 
lines represent the differential magnitude curves for the NB_1 .75 
and NB_2.12 filters, respectively. For comparison, the dashed 
curve shows the 5cr-sensitivity not applying LOCI but standard 
PSF subtraction to the 2.12 fim images. The LOCI algorithm 
improves the detection sensitivity by 0.5 up to 1.5 mag which 
corresponds to an improvement of up to factor of 4 in contrast 
ratio. The projected separation in AU was computed using the 
distance to the stars from Table [2] and is displayed at the bottom 
X-axis. The separation in arcsec is shown at the top x-axis. We 
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Fig. 3. 5cr-sensitivity profiles as a function of the distance from the star The NB_1.75 and NB_2.12 filters are shown in black and 
pink lines, respectively. For comparison, the dashed curve shows the 5cr-sensitivity for standard PSF-subtracted images at 2. 12 //m, 
not using LOCI. The cyan and green lines mark the approximate positions of 75 and 13.7 Mjup, the canonica l values that define 
the substellar and planetary regime boundaries according to the BT-DUSTY models from Al lard et ar| j 201CJ). T he dashed lines 
represent the same boundaries but for an age of 3 Myr for TW Hya recently reported bv lVacca & Sandel| (1201 ih . The gray areas 
show the disk inner hole/gap sizes reported in the literature and the red dashed vertical lines mark the intersections of the LOCI 
contrast curve at 2.12 //m with the cyan and green lines. They indicate the radial distance from the primary from where the NACO 
images are sensitive to substellar and planetary mass objects according to the adopted BT-D USTY model. The pink areas represent 
the hole/gap distance range where stellar companions could exist according to the models of lArtvmowicz & Lubowl (Il994h on basis 
of the hole/gap size. For typical gaseous viscous disk parameters, and a denoting the binary semi-major axis, the inner edge location 
of the circumbinary disk varies from 1.7a to 3.3a for binary eccentricity increasing from (circular orbit) to 0.75, and assuming 
a face-on geometry (/= 0°, PA= 0°). The position of the companion varies from a(l - e) and a(l + e), the periapsis and apoapsis 
of the secondary relative to the primary. For TCha, the NACO images are only sensitive to stellar mass objects brighter than 4.7 
mag in contrast or 0.12 M©. The blue vertical bar in the TCha plot indicate the approximate 5-cr detection limit {AKs > 4.65 mag 
or mass < 0.12 Mq) at the separation (6.7 ± I AU) and position angle (~ 78 + 1 deg) of the companion candidate detected in the 
L'-band with NACO bv lHuelamo et al.l(l20I ih . A lower limit of 50 Mjup was chosen following Fig. 4b in their paper This detection 
is well within the disk hole/gap of TCha and inside the binary range predicted by Artvmowicz & Lubow ( 1994). The profiles are 
displayed between 0. 1 and 1 arcsec for all stars. Simulations with a fake target of 4 mag in contrast show that the sensitivity curves 
are reliable as close as 100-120 mas from the primary (between the first and second Airy ring at 2.12 yum). Beyond 1 arcsec and up 
to 3.3 arcsec the contrast curves remain flat. 



typically reach a contrast of 2 to 5 mag at Of.' 1 and 4 to 6 mag at 
(y.'2. The mass limits can be estimated with a selected evolutio- 
nary model of a substellar object. The cyan and green lines mark 
the approx imate positions of the can onical values of 75 and 
13.7 Mjup (iBurrows et al.ll 19971 llooih that define the substellar 
and planetary regime boundarie^ for comparison with our data. 



The mass boundaries between stars, brown dwarfs, and planets are model-dependent. 
The former is defined as the minimum mass required for burning hydrogen for a solar metal- 
licity object in the main-sequence. This value is a function of the helium fraction, metallicity 
and opacity of grains and, depending on the assumed parameters, it can range from 0.07 to 



They were computed by interpolating the mass vs. NB_ 2.12 re- 
lation of the BT-DUSTY models from Allard et al. (2010f| con- 



0.092 Mo or 73 to 96 Mj„p IBurrows et al.ll997ll200il and references therein). The latter is 
defined as the mass-limit for deuterium burning and varies with the helium abundance, the 
initial deuterium abundance, the metallicity and the fraction of the initial deuterium abun- 
dance that must com bust for the substellar object to qualify as having burned deuterium 
ISpiegel et aU201Qh . The deuterium mass limit can range from ~ 1 1.0 to ~ 16.3 Mj^p de- 
pending on the model assumptions. A final remark is that other criteria, such as the formation 
process, should be considered when classifying a substellar object. 

These are the models of brown dwa rfs and ver y low-mass stars with dusty atmospheres 
fromlAllard et al. (2001) and Chabrier e t"ai] I2000l) but with updated opacities. One should 
be cautious when applying evolutionary models developed for stars and substellar objects 
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Table 2. Physical properties of the five stars and their disks found in the literature 



Star 


SpT 


K,' 


A/f,,' 


Dist. 
(pc) 


Ay 


Msiar 

(Mo) 


Age 
(Myr) 


(Moyr') 


DoAr21 


Kl" 


6.227 


0.018 


121.9±5.8'> 


6.2" 


1.8" 


<1" 


<io-"-' 


HD 135344B 


F4= 


5.843 


0.020 


140±42S 


0.5' 


1.7±0.2S 


8±48 


5xlO-''-s 


HR4796A 


AO' 


5.769 


0.023 


72.8±1.7'' 


- 


2.5' 


8±2' 


- 


TCha" 


GS'-'' 


6.954 


0.018 


108 ±9' 


1.7''-5 


1.5) 


7k,5 


4x10-" 


TWHya' 


K7" 


7.297 


0.024 


55±9'' 


0.01 


0.6' 


lOP, 3* 


8x10-'°-"- K)-'-' 


Disk 


Rdiiin-' 
SED, obs. 


Rhok/giip.in" 

SED, obs. 


Rhok/gap.oiir 

SED, obs. 


Rdi.k.o„,' 
SED, obs. 


/ 


PA 


Mui^k-* 


Class & disk type' 




(AU) 


(AU) 


(AU) 


(AU) 


(deg) 


(deg) 


(Mo) 




DoAr21 


- 


- 


-, 100" 


- 


<45'' 


- 


< 0.001' 


PMS circumbinary/debris' 


HD135344B 


0.18°. 0.05±0.25" 


0.45", 1.8±0.2"' 


45=, 39±4' 


300', 160S 


21'- 


55f 


(2.8±1.3)xlO-''8 


PMS Herbig e 


HR4796A 


- 


- 


40-200', 69' 


-,87' 


75.88±0.16' 


27.01±().I6' 


> 7.4x10-5-' 


Herbig debris disk' 


Telia's 


0.08«0.13",- 


0.2' 0.17°,- 


15' 7.5°,- 


300''',- 


75J, 60" 


78!^^-° 


0.003'°, (1.76±0.25)xlO-2.° 


WTT'^ 


TWHya' 


-0.02', - 


0.25', 0.06" 


4', 4.3±0.3" 


230" 


4.3±1.0''-',7±r-^ 


332±10'' 


-0.06' 


CTT* 



from 2MASS 

average accretion rate derived from Ha line einission for DoAr21 and TCha, Bry line emission for TWHya (y) and HD 135344B (g), and excess Balmer continuum for TWHya (q) 
radius from the primaiy derived from both SED fitting models and spatially resolved observations in the optical to the sub-mm/mm. In this work we adopt the latter values when 
available 

total disk masses estimated from mm fluxes assuming a gas-to-dust ratio of 100 or 77 for TWHya. TCha: the first value is that of the 3 mm emission flux for a distance of 66 pc. 
Because Mji^k ^ d", a value of 0.008 M© is expected for the assumed distance d=108 pc. The second value was derived from model fitting to the observed SED and H and A" -band 
interferometric observations with the VLTI/ AMBER instrument 

young stellar object (YSO) evolutionary class and disk type found in the literature: pre- main sequence (PMS ), weak lined T-Tauri (WTT) and classical T-Tauri (CTT) star 
TCha shows variable circumstellar extinction Av=1.2 - 4. 6 and an age of 4. 1 - 1 Myr ISchisano et aU2009l and references therein). Here we adopt the most frequent value Av= 1.7 
and age of 7 Myr. The disk inner hole size adopted is from lBrown et alj I2007D 

a recent work by '■ Vacca & Sandell 1 20 1 1 ) based on near-infrared spectroscopy revises the spectral type and age of TW Hya to M2.5 (instead of K7) and - 3 Myr (instead of 8-1 Myr). 
lAkeson et alj ^20 1 1 ) combines new near-lR inteiferometric data with previous spatially resolved observations at 10 fim and 7 mm to constrain disk models based on a flared disk 
structure. They find TW Hya fits a three-component model composed by an optically thin emission from - 0.02 AU to - 0.5 AU, an optically thick -0.1 AU wide ring that is followed 
by an opacity gap, and an outer optically thick disk starting at 3.8 - 4.5 AU. Here we adopt the values of 0.5 AU and 4.3 AU for the inner and outer hole/gap radius, respectively 



References. a) IJensen et afUIOOgh b)ILoinard et aUtlOOl. cjICiezaet aljt20Toll2003) d)[Barv et aUtlOOl. e)[Brown et aljt2007l) . fl IBrown et aUtlOMl). p) | Gradv et alj 120091) and 
references therein, h)'van Leeuwer^ ^2007^, i)'Schneider et al.' M999','200?').'Jura et al.' n99?),'Besla & WJ ^2007^ and references therein, j)'Alcala et al.' n993V 1 lk)'Schisano et^al.' 
|2009). 1) Torres etal. (2008). m) Lommen et al. (2007). n) Olofsson et al. (2011). o) Herbig (1978). p) Webb et al. (1999). q) Herczea et al. (2009). r) Calvet et al. (200^ 
s) Hughes et al. (.2007.). t) Thi et al. i20ia) u).Eisner et al. i2006) v.l) Pontoppidan et al. i2008.) for the inner disk, v,2) .0ietal. (.2004) for the outer disk, w) .Fedele et al. i200; 
x) lweinberger et alj 12553) , v) IEisneret'an 120101) . z) IAlcalaet al]|T99g) 



volved with the NACO filter at the distances and ages listed in 
Table|2l We then subtracted the 2MASS Ks absolute magnitudes 
of the stars (at the assumed distance an d corrected for extinction 
using the interstellar extinction law of iRieke & Lebofskvll 19851 
with Rv= 3.1) to obtain the difference in magnitude expected 
for a companion of 13.7 and 75 Mjup. The NB_ 2.12 filter was 
selected because a better contrast is expected at longer wave- 
lengths. The approximation of NB_2.12 mag ~ A's mag is a good 
first estimate (within < 0. 1 mag) that does not compromise the 
results if we account for all the errors involved. A younger age 
for TW Hya of 3 Myr (i nstead of 8-10 Myr) has been reported in 
IVacca & Sandelll ( 1201 lb . We computed the substellar and plane- 
tary mass boundaries of TW Hya for an age of 3 Myr and plotted 
it as cyan and green dashed lines for comparison. For this case, 
the NACO images are sensitive to substellar mass objects within 
the whole distance range probed (U.'! to 7'.'0 from the star). The 
gray area in the plots represents the size of the inner hole/gap de- 
termined from previous studies (Table|2]i. Three of the five disks 
in the sample have an optically thick inner ring of sub- to nearly 
2 AU in radii reported in the literature. Her e we adopt the inne r 
hole/gap radii of 1 8 AU for HP 1353 44B (iFedele et al.ll2008l) . 
0.2 AU for TCha dBrownet al."2007^d a6AU for TWHya 
from the recent results of [Akeson e t al.l (1201 ll) . Nevertheless, 
these values are not relevant for the analysis in this paper because 
the NACO observations are not sensitive to these inner radii. The 



to derive the mass of planets. The mass of a freshly formed planet only depends on the 
accretion history and accretion rate and a planetaiy model should be considered instead. 



vertical red dashed lines give an indication of the disk distance 
range where the images are sensitive to the substellar and plane- 
tary regimes, according to the BT-DUSTY model adopted here. 
The inner 100 mas of the profiles is omitted. The most inner 
central region of the residual maps is not dominated anymore by 
speckle or other noise but by systematic errors in removing the 
PSF core. In this region the "standard" contrast curve calcula- 
tion does not apply anymore. Simulations with a 4 mag contrast 
fake target show that the contrast curves are reliable to about 
the inner 100-120 mas or between the first and second Airy ring 
at 2. 12 jum. Table |2] lists the physical parameters of the five stars 
and their disks found in the literature. 



3.2. Orbit constraints on putative stellar/planetary 
companions 

This paper addresses one possible origin for inner holes or 
cleared gaps in protoplanetary disks, which is the dynamical 
clearing by unresolved companions. In this section, we discuss 
the expected orbits of potential companions that account for the 
observed hole/gap size on ba sis of existing mode l s or numerical 
simulations for both pla nets (Ouillen et al. 2004|; IVarniere et al.l 
l2006h and stars dArtvmowicz & Lubow..l994.) . 

We caution that these quantitative estimates are very uncer- 
tain because we do not know the orbital parameters and only 
poorly constrain the disk inclinations and other properties. The 
important result in this paper is qualitative and refers to the pre- 
sence or absence of stellar/substellar companions beyond the 0'.' 1 
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radius probed with the NACO images. We cannot completely 
rule out companions as long as there is a small chance of align- 
ment during the observations that could have occulted the com- 
panion. But the likelihood of that happening is small. Table 3 
shows the results for the distance range around each of the stars 
where we can effectively rule out the presence of stellar and sub- 
stellar companions on the basis of our NACO observations and 
with the limitation that we could have missed companions be- 
cause of unfavorable projections. 

Planetary companions 

A single planet of a few Mjup carves only a narrow gap (typi- 
cally Argap/rgap ~ 0.1) about its orbit in a viscous disk that is a 
few times more massive than the planet an d as soon as the Hill's 
radiu jj is larger than the disk height (e.g.. lLubow & D'Angeld 
I2OO6I) . This is not compatible with the large inner hole/gap 
observed in most transitional disks. Nevertheless, a system of 
multiple planets could open large holes and m aintain sharp in- 
ner disk edges (iPerez-Becker & Chiand 1201 iL and references 
therein). Another possibility is to consider that the mass of the 
single planet is large compared to the disk mass with which it 
interacts so that the inertia of the planet slows the migration 
and the inner disk within the planet's semi-major axis has time 
to viscously accrete onto the star, creating a large inner hole 
dOuillen et al.l 120041) . The material beyond the planet's orbital 
radius is prevented from accreting by the transfer of orbital an- 
gular momentum from the planet to the disk and piles up, form- 
ing the inner rim of the outer disk. Additionally, if the mas- 
sive planet forms faster than the migration timescale and slower 
tha n the hole/gap at its formation radii, numerical simulations 
by IVarniere et al.l (l2006h show that the "spectral hole" always 
accompanies the planet and extends from the planet's orbital ra- 
dius to the stellar surface. Here "spectral hole" subtends a region 
with reduced density, which is not necessarily fully depleted 
of material, allowing for accretion to hold even after the hole 
has formed. If the planet forms faster than the hole formation 
timescale, a gap can open with an inner radius larger than the star 
radius and an outer edge beyond the planet 's orbital radius. Other 
mechanisms such as outward migration (iMasset & Papaloizoul 
l2003h could contribute to a gap opening beyond the planet's or- 
bital radius. In this paper, and for simplicity, we assume that a 
planet of a few Mjup could exist within the whole extension of 
the hole/gap region of the five transitional disks in our sample. 



Stellar companions 

Binary systems with separations smaller than the outer disk, on 
the other hand, are expected to open inner holes that are signi- 
ficantly larger than the semi-major axis of the orbit of the binary, 
while keeping circumprimary disks that dissipate quickly on vis- 
cous timescales. In this case, the inner hole or gap is created by 
tidal truncation of the disk which, initi ally interpreted as tran- 
sitiona l, is in fact a circumbinary disk. lArtvmowicz & Lubowl 
(1 1 9941) investigated the gravitational interaction of an eccentric 
binary with the circumbinary and circumstellar gaseous disks 
and analytically computed the approximate size of the disk 
gaps as a function of the binary mass ratio and eccentricity. 
For typical viscous disk parameters and a denoting the semi- 
major axis of the binary, the inner edge location of the cir- 



cumbinary disk relative to the center of mass of the binary, or 
Rhoie,<mi, varies from l.la to 3.3fl for a binary eccentricity in- 
creasing from (circular orbit) to 0.75 and assuming a face-on 
geometry (/ = 0°, PA = 0°). Therefore, R hole, out I ^■'i < a < 
R hole, out I ^■'^ and the position of the companion ds varies from 
a{\ - e) to fl(l + e), the periapsis and apoapsis of the secondary 
relative to the primary. Observations of binary systems inside 
disk holes suggest that these numbers are approximately correct 
(e.g., iGuenther et aP IIoOTI for CS Ch a. Ilreland & Krausl I2OO8I 
for Co kuTau/4.lBeust & Diitrevll2006l for GG Tau, and more re- 
cently, |Huli^oetal]|20Tl] for TCha). The circumbinary disk 
recedes from the binary with increasing eccentricity. For sim- 
plicity, we assume the origin of the coordinate system to be co- 
incident with the primary star This will not affect the results for 
ds significantly because the clearing is always larger than a and 
because of the many uncertainties involved. 

In order to compute the expected positions of stellar com- 
panio ns for the stars in our sample , we considered the exam- 
ple in lArtvmowicz & Lubowl d 19941) of a binary star with vis- 
cosity parameter a^ ~ 10"^, Reynolds number R ~ 10^ and 
binary mass parameter// = 0.3 for which Rhoie.om = l-7fl for 
e = (model I), Rhoie,oui - 2.6fl for e - 0.25 (model 2), and 
Rhoie,oui - 3.3a for e — 0.75 (model 3). The results are summa- 
rized in Table 3 and are plotted in Fig. 3 as the pink area in the 
sensitive curves. This area represents the extension within the 
disk hole/gap where stellar companions could exist according to 
these models and for a face-on geometry. The upper limit of this 
region is given by the circular orbit (model \,ds - Rhoie.oui/^-^) 
and the lower limit by the periapsis of the most eccentric orbit 
(model 3; ds = fl(l - e) = Ri,oie_„,„/ 3. 3(1 - 0.75)). The cor- 
responding orbital period P was computed assuming Keplerian 
rotatiorQ and m - 0.5M* at the location where companions are 
expected. 



Projection effects 

In this section we evaluate the orbit projection effects on our 
results on the basis of the disk inclinations / reported in the lite- 
rature for the five sources, assuming the binary is coplanar with 
the disk. Companions in nearly edge-on systems are likely to be 
occulted by the primary or by the disk itself. The gap properties 
also remain unchanged unti l the disk is nearly edge-on and starts 
to obscure the central star ( iBrown et alj|2007l) . Assuming a cir- 
cular orbit for simplicity (we could not access the eccentricity 
of a detected companion with just one epoch images, anyway) 
and polar coordinates with the origin of the referential at the pri- 



mary, the projected orbit can be computed with Xproj = a cos 6 
and yproj - a sin 6 cos /, where a is the semi-major axis of the 
real orbit, / the inclination and the angle between the position 
vector and the x-axis. The projected distance between the pri- 
mary and secondary is determined by dproj = ix~ . ■ + y^ . .)'^~ 
at any point in the orbit. The results are listed in Table 3 and, as 
expected, the projection effects are more significant for the more 
inclined disks HR4796A and TCha. 

In our case, a non-detection of the companion would occur 
for projected distances shorter than (X'l, the detection limit in 
radii of the NACO images. Confining the analysis to the first 
quadrant (0° < < 90°), we can determine the angle for which 
dproj < 0'.' 1 and predict the percentage of the orbit we could ef- 
fectively detect with our NACO observations. We performed this 
exercise for two values of the semi-major axis, a - Rhoie.om / 1 -T , 



The Hill radius of a planet is given by i?// = ap{l — e){j^y^^, where ap is the semi- 
major axis of the orbit of the planet, e the eccentricity, nip the mass of the planet and M* the 
mass of the star. 



and m the mass of the secondary. 



where G is the gravitational constant. A/* the mass of the primary 
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Table 3. Predicted binary signatures from the models of lArtvmowicz & Lubowl (Il994t) . projection effects caused by the incUnation 
of the orbit, and hole/gap range where stellar/substellar companions can be ruled out from the NACO images 



Star 


RholL-,oul 


Model I' 


Model 2- 


Model 3-' 


/ 


Pi-oj. etfects-* 


ARsldlar^ 


AR.ubs,cii>„-' 




ds 


P 


ds 


P 


ds P 


dp,„j* dp,„i" 






(AU) 


(AU) 


(yr) 


(AU) 


(yr) 


(AU) (yr) 


(deg) 


(AU) (AU) 


(AU) 


(AU) 


DoAr21 


100 


59 


276 


29-48 


95 - 202 


7.6-53 13-235 


<45 


42-59 71-100 


12.2-100 


31 - 100 


HD135344B 


39 


23 


69 


11-19 


23 - 52 


3-21 3-60 


21 


21-23 36-39 


14-39 


37-39 


HR4796A 


69 


41 


136 


20 - 33 


46-98 


5-37 6-116 


76 


10-41 17-69 


7.3 - 69 


13-69 


TCha 


15 


8.8 


17.4 


4.3-7.2 


6-13 


1.1-7.9 0.8-15 


60 01- 75 


2.3-8.8 4-15 


10.8- 15 


- 


TWHya 


4.3 


2.53 


4.2 


1.25-2.1 


1.5-3.2 


0.33 - 2.3 0.2 - 3.7 


4.3 


2.52-2.53 4.29-4.3 


- 


- 



Model 1 (circular orbit): e — 0, Rhote.out - 1 .7fl and the distance of the secondary relative to the primary is ds 
Model 2: e = 0.25, Ri,„ir,„,„ = 2.6a and a(l - e) < ds < fl(l + «) or 0.29R;,„fe,,„„ < ds < 0.48R;,„/,.„„, 
Model 3: e = 0.75, R,,ah.„„i = 3.3fl and a(l - e) < ds < a(\ + e) or 0.076«,,„fc„„, < ds < 0.53«,,„fc„„, 



^0), 



: a sin Gcosi and d,„ 



„n.i 



- (■'^Iroi +>'pra/)"''''p'"j(""'-'") = "for" = 0° andfl',,„,j(Him) = a cos i 



Binaiy coplanar with the disk {is = i^i.^kX circular orbit (e - 

for = 90° . Computed for *« - Rhcie^ml 1 -"^ and **fl - Rijoie^m 

Range of disk hole radii, starting at (K' 1 from the primary, in which stellar binary companions can be ruled out on basis of the VLT/NACO images and using the BT-DUSTY 

evoludonary models 

Range of disk hole radii in which suhslellar binary companions can be ruled out from this work 



and a - Rhoie,our- The results show that a non-detection of the 
companion owing to projected distances shorter than (J/ 1 could 
exist only for TCha for a = Rhoie.out = 15 AU. At a distance of 
108 pc ((y.'l ~ 10.8 AU) and for disk inclinations / = 75° or 60° 
reported in the literature, a non-detection could occur if = 45° 
or 50° (1st quadrant) or when 45 - 50° < < 225 - 230° and 
225 - 230° < 6* < 315 - 320°, for the entire orbit. Therefore, 
we would access 50 to 55% or 19.5 to 21.5 yr of a total orbit of 
39 yr for a distance of 15 AU from the primary and m - 0.5M*. 
If a = Ri,oic/l.7 - 8.8 AU, the NACO images are no longer sen- 
sitive to the distances involved. This is also the case for TW Hya, 
for which a.'l ~ 5.5 AU and Rnoh.oui = 4.3 AU. 

4. Discussion: hole/gap formation mechanisms 

The origin of transitional disks is still under debate with seve- 
ral theories developed to explain inner opacity holes and gaps in 
protoplan etary disk s, including _i) magn eto-rotational instabili- 
ty (MRI: Icidang & Murrav-ClayJ 
by the central s tar dClarke et all U 



2007h. ii) photoevapo ration 
20011: [Alexander et alJl2006t 



Gorti et al.l2009l). iii) dust rernoval b y coagulation jTanaka et alj 
2005t iDuUemond & Dominikll2005h and iv) dynamical clearing 
by un resolved companions , both stars ([Artvmowicz & Lubow 
1994), and giant planets dOuillen et al.1 l2004t IVarniere et al. 
2004 . 

The MRI mechanism operates in ionized gas and hence, at 
all disk surface layers ionized by direct X-rays from the cen- 
tral star, including the inner disk wall. Although this still al- 
lows accretion onto the star, MRI cannot explain the remnant 
optically thick disk observed in "pre-transitional" disks. Current 
photoevaporation models require a maximum total disk mass of 
~ 5 Mjup and negligible accretion rate, unless the disk is un- 
usually inviscid or large (> 100 AU). Extreme-UV (EUV), far- 
UV (FUV) photons and X-rays from the central star heat and 
ionize the circumstellar hydrogen, generating a disk evapora- 
tive wind beyond some critical radius. At disk early-evolutionary 
stages, the accretion rate dominates over the evaporation rate. 
For photoevaporation to hold, the accretion rate has to fall below 
the wind rate and, once this happens, the photoevaporative wind 
can open a gap in the gas disk. The dust in the inner disk is then 
rapidly removed (on a viscous timescale of typically less than 
10^ yr), forming the inner hole, and the entire disk dissipates 
very quickly through photoevaporation. Computed disk lifetimes 



are ~ 4 M yr for M» g 3 M p and less than 1 Myr for more mas- 
sive stars ( iGorti et alj2009t) . The inner disk drains away entirely 
before the gap can grow beyond l-IO AU, and therefore photo- 
evaporation cannot explain the inner dust disks and large inner 
holes observed in many transitional objects. Grain growth can 
also produce opacity holes but is a strong function of the ra- 
dius. It is more efficient at disk small radii where the surface 
density is higher and the dynamical timescales are shorter and 
do not affect accretion (of the gas). Unfortunately, is not possi- 
ble to detect grains larger than a few millimeter and, therefore, 
this mechanism was not observationally confirmed in the last 
years. And finally, as discussed in section \T7\ stellar/substellar 
companions or giant planets of 1 to a few Mjup can open gaps 
and inner holes through tidal interaction with the disk. While a 
stellar companion is expected to almost suppress accretion, an 
embedded planet of ~ 1 Mjup has only a modest effect on the 
disk accretion rate, reducing it by a factor of < 10 inside the 
planet's orbit compared to that from outside, and ha s little effect 
on the total disk mass dLubow & D'Angeloll2006h . And, for a 
given disk mass, the more massive the pl anet and the less viscous 
the d i sk, the lower the ac cretion rate ([Alexander & Armitagd 
l2007h . iLubow et all (Il999h proposed that a planet of 10 Mjup or 
more could halt stellar accretion. Nevertheless, a few observed 
close stellar/s ubstellar binaries have been reported to be accret- 
ing: DQTau dCarr et alJ l200ih. CS Cha (lEspaillat et alJ l2007h 
and, m ore recently, T Cha dHuelamo et al.ll201ll : lOlofsson et al.l 
1201 ll) . Factors other than the presence and mass of a close com- 
panion must determine accretion onto the star and hence, the ac- 
cretion rate alone cannot be used to infer or discard the presence 
of potential companions. The viscosity and scale height of the 
disk and eccentricity of the orbit may contribute to maintain the 
accretion across th e disk gap/hole even for stellar companions 
d Artvmowicz & Lubow 1996) . 

[Alexander & Armitagd d2007h proposed a simple observa- 
tional diagnostic to distinguish between photoevaporation and 
giant planet-induced inner holes that relies on the relative values 
of the the total disk mass (Mdisk) a nd stellar accretion rate (Mace). 
More recentlv. lCieza et aP d2010l) used optical-to-millimeter ob- 
servations of a sample of 26 transitional disks in the Ophiuchus 
molecular cloud (d ~ 125 pc) to derive the disk parameters 
and investigate the mechanisms that potentially create their in- 
ner holes. In addition to the disk mass and accretion rate, the 
authors also consider the fractional disk luminosity (Ldisk/^*) 
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and parameters from the SED. According to their criteria, 
grain-growth-dominated disks have high-accretion rates (Mace ^ 
10"** Moyr"') and large disk masses because this mechanism 
mostly operates in the disk inner region, whereas most of the 
mass resides at larger radii. Photoevaporating or photoevapo- 
rated disks should have a small disk mass (Mjisk < 2.5 Mjup) 
and negligible accretion (< 10 " Moyr '), and if in addition 
the fractional disk luminosity Ljisk/i'* < 10"^, the disk could 
be already in the debris evolutionary stage regardless the ori- 
gin of the inner hole. Good candidates for ongoing giant planet 
formation are disks with large mass (1.5 - 11 Mjup) and low ac- 
cretion rate (~ 10 '" - 10"^ Moyr"'), while disks with small 
mass (< 1.5Mjup) and higher accretion (lO"'^-' - lO"^-' Moyr"') 
could have formed a giant planet in a recent past because most 
of the disk mass has alread y been depleted. Cl ose binaries are 
ruled out on the basis of the lArtvmowicz & Lubow (il994i) mo- 
dels (Rhole.oul ~ 2a). 

Following these criteria and using the disk parameters in 
Table 2, we tentatively classify the five transitional disks in our 
sample on the basis of the mechanisms that potentially create 
their inner opacity hole/gap. DoAr21 was classified as a pos- 
sible circumbinary or debris disk by .Cieza_et al,: (20T3) (source 
12 in their catalog). HD 135344B, TCha, and TWHya are good 
candidates for ongoing giant planet-formation because of their 
large disk masses and modest accretion rates. HR4796A could 
be a photoevaporated, debris, or a planetary disk owing to its 
very low mass, dust ring structure and no reported accretion rate. 
Nevertheless, this classification is just a best guess given the 
available data and the uncertainty associated to all parameters, 
both observational or derived from the models. Additionally, 
the hole/gap clearing mechanisms are not all mutually exclu- 
sive. Two or more mechanisms can operate at the same time 
(for example, grain growth and MRI, or planet formation and 
photoevaporation) or sequentially (once the gap has formed by 
one mechanism, others can rapidly remove the inner disk ), con- 
tributing to the present observed transitional evolutionary state 
of the disk. Transitional disks seem to represent a highly hete- 
rogeneous group of objects exhibiting a wide range of masses, 
accretion rates, and SED morphologies although, overall, they 
tend to have s maller masses and lower accretion rates than disks 
without holes (ICieza et al.ll2010l) . To conclude, disk dissipation 
does not seem to follow a standard evolutionary path, but instead 
many different paths are possible to which different mechanisms 
contribute. 

The VLT/NACO images presented here provide an observa- 
tional input to the last proposed mechanism by ruling out the 
presence of stellar companions within part (beyond O'.'l) of the 
inner hole/gap region of four disks in our sample and down to 
the planetary mass regime for DoAr21 and to the substellar 
mass regime for HD 135344B and HR4796A. The exception is 
TW Hya, for which the NACO images cannot resolve the inner 
hole of just 4.3 AU, or 78 mas for a distance of 55 pc. The re- 
sults in Fig. 3 (pink bars in the plots) and Table 3 show that our 
observations are sensitive to part of the region where binaries 
are exp ected to exist from the models of lArtvmowi cz & Lubowl 
d 19941) for HD 135344B and almost this entire region in the case 
of HR 4796 A and DoAr21. With the proviso of possible occul- 
tations, we do not report any detection of companions beyond 
O'.'l from the primary for the five transitional disks. Therefore, 
the results in this paper further constrain the suite of potential 
explanations for their origin and toward the exciting possibili- 
ty that giant planets are currently forming inside their inner 
holes/gaps. Nevertheless, the fraction of transitional disks that 



are in fact tight binaries still remains to be established. Common 
adaptive optics observations (on 8 to 10m class telescopes) can- 
not constrain the inner 0'.' 1 region from the primary, or radii of 
less than 5 to 15 AU for most transitional disks. Additional sur- 
veys using different techniques such as sparse aperture masking 
(SAM) with adaptive optics instruments (the one leading to the 
detection of the potential substellar companion of T Cha reported 
in iHuelamo e t al. 201 1) and radial velocity studies, should be 
performed to probe the inner most regions and firmly establish 
the fraction of transitional disks that are in fact close binaries. 
Finally, the Atacama Large Millimeter Array (ALMA) with its 
1-3 AU spatial resolution and high sensitivity will be capable 
of resolving young planetary systems, detecting and determi- 
ning the properties of close companions (both planets and stars) 
and of the dust and gas through out the disk. Planet and inner 
hole/gap formation scenarios will be confirmed or ruled out. 

The following sections summarize the results regarding the 
presence of potential companions for each individual source 
from previous studies and this work. 

DoAr21 



iJensen et al.l (l2009h summarizes the properties of this object. 
Gemini sub-arcsecond resolution 9-18 yum images show little 
or no excess mid-IR emission within 100 AU from the star. 
But beyond that the emission extends over several arcsec (or 
hundreds of AU) and is quite asymmetric, showing a half -ring 
structure of higher emission at I'.'l from the star to the NW. 
This feature is bright in PAHs emission at 8.6 and 11.3 //m, 
H2 emission at 2.12 jjm (the 1-0 S(l) ro-vibrational line) and 
continuum emission at 18 fim, and very different from what is 
expected from a disk. DoAr21 is a particular object different 
from other pre-main sequence stars (PMS) because it is very 
young, a very strong emitter of (hard) X-rays and shows no ac- 
cretion. Based on the estimated FUV and X-ray flux at the po- 
sition of the ring (~130 AU) and the prop e r moti on of DoAr21 
(9.5 + 1.3 kms"' westward). iJensen et al.l (l2009l) proposed this 
emission to be associated with a small photodissociation region 
(PDR) resulting from the interaction of the radiation field with 
denser material from the Ophiuchus cloud in the vicinities of 
the star and not with a circumstellar (or circumbinary) disk. The 
half-arc shape would form because the material is swept up by 
the stellar motion. The H2 emission was spatially resolved with 
the AO-assisted IFU instrument SINFONI at the VLT as an arc 
of material circling the star at radii between 0'.'6 (73 AU) and r.'8 
(220 AU) and position angles of -40° to 170° by M. Hogerheijde 
et al. (to be submitted), who proposed several explanations for 
the asymmetric emission. The derived H2 line ratios were con- 
sistent with thermally excited gas by X-rays (T = 1000-2000 K) 
and not with FUV fluorescence. The most likely scenarios in- 
clude illumination of unrelated cloud material or material resul- 
tant from the collision between two (proto)planets and the mi- 
gration of a substellar/planetary companion to the ring region 
(73-220 AU) that could disrupt a previously stable disk. The last 
scenario was ruled out by their SINFONI images down to a mass 
limit of 5 Mjup, by assuming an age of 0.3 Myr for DoAr21. 
However, an unseen planet of a few Mjup could still explain the 
observations. Future epoch observations are necessary to deter- 
mine if this material is orbiting the star (in a disk) or if is tra- 
veling through the cloud together with DoAr21, confirming the 

PDR scenario. 

VLBA observations bv lLoinard et al.l (l2008h show DoAr21 
as a binary with 5 mas (0.6 AU) projected separation and esti- 
mated semi-major axis of 1-2 AU, but no information is given on 
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the mass or optical luminosity ratio. Assuming it to be an equal- 
mass binary, evolutionary tracks predict for each star a mass of 
~ 1.8 Mq (~ 3.6 M0 for the binary) and an age of ~ 8 x 10^ yr. 
Nevertheless, such a close companion is beyond the resolution 
ofNACO. 

As expected, the H2 emission was not recovered in our 
VLT/NACO 2.12 /dm images (the width of the H2 Hne in this 
source is of 2x10""* yum and was detected with an integration 
time of 7 seconds and a S/N of ~10). Taking into account the 
width of 0.022 jum of the NB^.12 NACO filter, one would 
have needed integration times 220 times longer (i.e., 25 min) 
to achieve a similar S/N with our observations. Despite failing 
to detect off'-source emission at 2.12 yum, the NACO images 
confirm the absence of substellar/planetary companions in the 
arc region down to a mass of 15 Mjup at 73 AU and 6 Mjup at 
192 AU. Our NACO images conservatively exclude substellar 
objects (m < 75 Mjup) at radii larger than ~ 31 AU (0'.'25) and 
of stellar companions within the extension of the circumstellar 
environment probed with these observations (r > CK'l or 12 AU). 
The 5cr sensitivity at (X'l is 3.1 mag or 0.19 M0 for an assumed 
ag e of 1 Myr. 

ICieza et aP (l2010l) collected nearly diffraction 
VLT/NACO/S13 images in the J- and iT.-bands (FWHM 
of 60-70 mas) from April to September 2009 of several 
transitional disks, including DoAr21. A visual inspection of 
those images revealed no stellar companions beyond ~ 8 AU 
from the primary, which agrees with our results. Their detec- 
tion limits at (X'l separation, estimated from the 5cr noise in 
the PSF subtracted images, have a median value of 3.1 mag 
corresponding to a flux ratio of 17. These authors classified 
DoAr21 as a possible circumbinary or debris disk based on 
the disk parameters and no observed accretion ( see discussion 
above). As proposed also bv Ijensen et al.l (l2009l) . the DoAr21 
disk must have almost entirely dissipated either through binary 
formation or photoevaporation. This scenario is indeed possible. 
An equal-mass eccentric binary with semi-major axis of 1-2 AU 
would tidally clear out an inner hole of 3 to 6 AU. Once the gap 
was opened, the high FUV and X-ray luminosity of the binary 
with mass of ~ 3.6 M© w ould photoevapor ate the disk on a 
timescale of ~ 10^ - lO*" yr ( iGorti et alj|2"009l) . 

The NACO images of DoAr21 are the only ones from the 
sample of five disks that are sensitive enough for a detection of 
objects less massive than ~ 13 Mjup that is, potential planets or 
low-mass brown dwarfs within the inner hole region of the disk 
and at radii larger than ~ 76 AU (0'.'63). 



HD 135344B or SAO 206462 

Di rect imaging a t 880 ;um with the Submillimeter Array (SMA) 
by iBrown et al.l (l2009h spatially resolved the disk and an inner 
hole extending to a radius of 39 AU, for a distance of 140 pc 
to HD 135344B, in agreement with the value of 45 A U de- 
rived from the SED best-fit model (Brown etal. 20071 ). Thi s 
star has a wide companion at 20'.'4 CCoulson & Walthenll995h . 
HD 135344B has been observed by the HST in coronographic 
mode with NICMOS (1.1 and 1 .6 yum) and STIS and no com- 
panions were found dGradv et a l. 2005, 2009). The 5cr detec- 
tion sensitivity limits show the non-existence of stellar compa- 
nions within ~ 0.1 to 1 a rcsec (~ 14 to 140 AU) from the star 
iPontoppidan et al.l (|2008|) detected CO emission at 4.7 ^m ex- 
tending from 0.3 to 15 AU from the primary with CRIRES at the 
VLT, arguing against stellar clearing, but in favor of a planetary 
mass companion o rbiting at 10-20 AU. Spectroastrometry by 
iBaines et al.l (l2006l) also excludes stellar companions farther out 



than 0.42 AU. All these studies cannot exclude substellar objects 
within the disk (160 AU). The current NACO observations ex- 
tend the previous results by excluding brown-dwarf-mass com- 
panions at orbits larger than ~ 37 AU (0'.'27) and down to a mass 
Umit of ~ 19 Mjup (8 mag in contrast) at 221 AU. At O'.'l they 
are sensitive to stars more massive than 0.22 Mq. 



HR4796A 



[Schneider et al.l ( |2009|) obtained high spatial resolution 
(~ 70 mas) optical coronographic images of the HR4796A 
debris dust ring with the Space Telescope Imaging Spectrograph 
(STIS) resolving emission from -69-87 AU (0'.'95 - 1'.'2) in 
radii (see Fig. 7 in their paper) and a displacement of 1.4 AU 
between the disk center and the star, measured along the ring 
major axis. These values are in accordance with previous obser- 
yations at 1.1 and 1.6 /urn with the HST/NICMOS instrument 
dSchneider et al.lll999l). The hole ou ter radius was first inferred 
from the SED bv iJura et al.1 (1 19951) to be from 40 to 200 AU. 
HR4796A is a visu al binary having a M2.5-dwarf companion 
at 7'.'7 or 560 AU (iJura et al.lll993l) . More recendy, a tertiary 
wide companion (~ 3'separation) was disc overed from its X-ra y 
emission and is a M4.5 spectral-type star dKastner et alJl2008h . 
No other companions have been detected in H- and ^j-band 
coronographic NA CO imaging at proj ected separations in the 
range CK'l to 10" dChauvin et al.ll2010l) an d there are n o radia l 
velocity measurements for this source. IWvatt et all d 1 9991) 
concluded that a planet more massive than 10 MEaith close to 
the inner edge of the ring would be required to explain the 
asymmetry seen in the dust emission. Such a low-mass object 
could not be detected with the observations presented here. 
Nevertheless, the NACO images discard the existence of stellar 
companions within the disk hole at radii larger than ~ 7 AU 
((X'l) and more massive than 0.24 Mq, and of substellar objects 
outwards of ~ 13 AU (O'.'l 8). The mass detection limit at the 
hole outer radius (69 AU) and inside the dust ring (up to 87 AU) 
is 17 Mjup. 



TCha 

Radial velocity (RV) observations we re performed f or thi s 
source and the results are presented in ISchisano et al.l d2009l) . 
The authors report a radial velocity variation of 10 kms"', albeit 
non-periodic. No companions for TCha were detected at pro- 
jected separations of (X'l to 10" in diffraction-l imited H- and 
ATj-ba nd NACO images wit h Lyot-coro nography dChauvin et al.l 
l2010b . Recentlv. lHuelamo et al.l d201 lb reported the presence of 
a potential substellar companion at a separation of 62 + 7 mas 
(6.7 + 1 AU) from TCha detected in the L'-band (contrast AL' - 
5.1 + 0.2 mag) with NACO and using sparse aperture masking 
(SAM). The object was not detected in the authors' /iTj-band data 
up to a 3cr contrast limit of 5.2 mag. The corresponding 5cr con- 
trast Hmit is ~ 4.65 mag or 0.12 Mq for the BT-DUSTY model 
adopted here and is indicated in Fig. 3 as the blue vertical bar in 
the T Cha sensitiv ity plot. A lower limit of 50 Mjup was chosen 
following Fig. 4 in JHuelamo et al.ld201 ll) paper. This detection is 
well within the TCha disk hole/gap and inside the binary range 
predicted by the models of lArtvmowicz & Lubowl dl994|] i 



Assuming the companion is coplanar with the disk, which has an inclination of 60" to 
75° repotted in the literature, and a circular orbit, we can put some constraints on the orbit 
of the potential companion by considering the two extreme cases: 1)^ — 0° and the distance 
primary-secondary equals the binary semi-major axis, d^ = a = 6.7 AU, and 2) 6 = 90° and 
J,i = acosi = 6.7 AU. In 1) the potential companion is at its largest projected separation 
at the moment and will reach the shoitest separation of 1.7 to 3.3 AU when 9 = 90° for 
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Although we could not have detected such an i nner low-mass 
object, these observations support ours and the IChauvin et al.l 
(1201 Ol) results by ruling out stellar companions outward of (X' 1 
from the primary. The NACO images of T Cha exclude stellar 
companions brighter than 2.2 mag in contrast (0.29 Mq) at a 
projected separation of O'.'l to 2.9 mag (0.25 Mq) at 15 AU, or 
within the disk inner hole. They are not sensitive enough for the 
detection of substellar binaries within the whole extension of the 
disk. The maximum attainable contrast is 4.7 mag, correspon- 
ding to a 5cr sensitivity of 0.12 M©. 

TWHya 

Detailed modeling of the SE P of the star pred icted a disk 
hole/gap outer radius of ~ 4 AU (ICalvet et al.l2002h later on con- 
firmed with Very La rge Array (VLA) observations at 7 mm by 
iHughes et alJ (l2007h . The data resolved the hole in the dust disk 
and revealed a bright inner rim asso ciated with the disk e dge 
being directly illuminated by the star. ISetiawan et al.l (l2008h re- 
ported a 10 Mjup companion at 0.04 AU from the star through RV 
observations in the optical regime. jHuelamo et a l. ( 2008|) car- 
ried out near-IR high-resolution observations of the source with 
the CRyogenic high-resolution InfraRed Echelle Spectrograph 
(CRJRES) mounted on the VLT and found no trace of this com- 
panion. Spectro- astrometric imaging of the 4.7 fixn rovibrational 
lines of CO by iPontoppidan et al.l (12008!) detected gas in the 
TW Hya disk hole between 0. 1 and 1 .5 AU. Based on this re- 
sult, the authors suggest that the presence of a ~16 Mjup body 
could be clearing the dust and rule out photoevaporation as the 
hole formation mechanism. They further explain that if the 4 AU 
hole is created by a planet, it should be placed in an orbit of 
1-2 AU. The present NACO observations cannot resolve these 
inner disk regions. 

Direct coronographic imaging of the disk in scattered light 
at 1.1 an d 1.6 urn with the HST/N ICMOS instrument was col- 
lected by [Weinberger et al.l (l2002h . The scattering profile indi- 
cates a flared disk extending from ~20 to 230 AU. No point 
sources were detected with a 3cr significance from ~ 0'.'54 to 
4" (or 30 to 220 AU) down to a mass of ~ 10 Mj^p at ~ 0'.'63 
(35 AU) and 3 Mj^p at 1'.'52 (85 AU). Additionally, TW Hya was 
included in an a daptive optics imaging s urvey with the Keck II 
10m telescope bv lBrandeker et all (l2003h . The star was observed 
in the //-band and reached a 5cr contrast sensitivity of ~ 2 mag 
at a separation of 0'.' 1 . No companions were detected between 
^.'042 (the diff'raction Hmit in the //-band) and 1'.'6, or from 2.3 
to 88 AU. 

The results derived here from the NACO imaging further 
complement and support the previous studies by excluding stel- 
lar binaries more massive than 0.11 M© (3 mag in contrast) out- 
ward of ~ 5.5 AU from the primary (the inner detection limit 
of the observations or 0'.' 1), brown-dwarf -mass companions out- 
ward of 7 AU (0*.' 13) and planetary mass companions outward of 
~ 13 AU (0'.'24). We reach the maximum sensitivity of 7 Mjup, 
or 8 mag in contrast fr om 87 AU onward . 

A recent paper by IVacca & Sandelll (1201 ih based on near- 
infrared spectroscopic data revises the spectral type and age of 
TWHya to M2.5 (instead of K7) and ~ 3 Myr (instead of 8 - 
10 Myr). Under the assumption of a younger age, the NACO 



an inclination varying between 60° and 75°. The orbital period varies between 13.6 yr (for 
m = 0.12 Mq) and 14 yr (for in = 50 Mj^p). In 2) the companion is at its shortest projected 
separation in the present and the semi-major axis of the orbit a can be derived assuming an 
inclination of the orbit /: for / = 60°, a = 13.4 AU and for / = 75°, a — 26 AU. The orbital 
period varies between ~ 39 yr (/ - 60°) and ~ 105 yr (/ - 75°). Future epoch observations 
will better constrain the orbit of the potential companion by providing the inclination and 
the eccentricity. 



images are now sensitive to substellar mass companions within 
the whole disk region probed (0'.' 1 to 7'.'0) and to planetary-mass 
objects beyond ~ 10 AU (a.' 18). 



5. Conclusions 

This paper addresses one possible origin for inner holes or 
cleared gaps in protoplanetary disks. These holes/gaps are seen 
as potential signposts of planet formation but other explanations 
are possible. Here we investigate the mechanism of tidal trun- 
cation by close binaries as the possible or alternate explanation 
for the origin of the inner holes/gaps observed in five transitional 
disks. The conclusions are that: 

1. With the proviso of possible unfavorable projections, the 
VLT/NACO images rule out the presence of binary compa- 
nions from O'.'l to ~ 7" in all five transitional disks and down 
to a 5cr detection limit of 2.2-4.9 mag in contrast for the in- 
ner region (at O'.'l) and 7.7-9.3 mag (4.7 mag in the case of 
TCha) for the outer region (beyond I'.'O). 

2. For the transitional disks DoAr21, HD 135344B, HR4796A, 
and T Cha, the VLT/NACO observations resolve part of the 
inner hole/gap region of the disk and can potentially discard 
binaries down to a certain mass limit as the main explanation 
for the inner opacity cavities. For these objects, our results 
reduce the suite of potential mechanisms that create the holes 
to closer or lower-mass companions, giant planet formation, 
efficient grain growth or photoevaporation (for DoAr21 and 
HR4796A). 
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